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ABSTRACT 

The analysis of ion data from retarding potential analyzers (RPA’s) 
is generally done under the planar approximation, which assumes that 
the grid transparency is constant with angle of incidence and that all ions 
reaching the plane of the collectors are collected. These approximations 
are not valid for situations in which the ion thermal velocity is compar- 
able to the vehicle velocity, causing ions to enter the RPA with high 
average transverse velocity. To investigate these effects, we calculate 
current-voltage curves for H + at 4000° K, taking into account the finite 
collector size and the variation of grid transparency with angle. These 
curves are then analyzed under the planar approximation. The results 
show that only small errors in temperature and density are introduced 
for an RPA with typical dimensions; and that even when the density error 
is substantial for non-typical dimensions, the temperature error re- 
mains minimal. This last fact may be useful for designing a non-typical 
RPA which is optimized to operate in unusual conditions. 
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One of the primary instruments for making direct measurements on charged 
particles in the earth's ionosphere has been the retarding potential analyzer 
(RPA). RPA's have been flown on rockets, satellites and space probes to de- 
termine electron density and temperature, and positive ion density, temperature 
and mass number (Serbu and Maier, 3966; Maier, 1969; Chandra et al., 1970). 

A typical planar RPA is represented in Figure 1, with (for example) three 
grids G, , G 2 , and G 3 , and a current collector C. The aperture radius is r, the 
collector radius is R, and the interelement spaclngs are dj , d 2 , and d 3 . Charged 
particles entering the aperture and reaching C are collected as a measured 
current I, For ion measurements, constant voltages applied to G t , G 3 , and C 
repel the electrons and accept the ions, so that I consists of only ions. A varying 
ion retarding potential V is applied to G 2 , so that I is a function of V. Further 
details of RPA operation may be found elsewhere (Knudsen, 1966; Donley, 1969). 

The parameters to be measured are then obtained by fitting a theoretical 
curve to a large number of experimental data points comprising the collected 
current versus retarding voltage curve, or I-V curve. Tills fitting process is 
generally more involved for ions than for electrons, since more than one ionic 
species may be present, and the resulting I-V curve is the sum of the individual 
I-V curves for the several species. Despite this complication, programs have 
been developed which can deduce the ion temperature and the individual densities 
from either simple one-component curves or from multi -component curves 
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(Mobs and Hyman, 1968; Patterson, 1869), In these analysis programs, the 
I-V curve is assumed to obey the following formula, Witten here for multiple 
ion species : 
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= satellite velocity 
= trap-velocity angle (Fig, 1) 


= thermal velocity of i tl1 ion species 
= Boltzmann's constant 
= ion temperature 
« mass of i t!l ion species 
= electron charge 
= retarding potential 
= transparency of single grid 
= ion density of i tl1 ion species 
= satellite potential. 



(i) 


All ion species are assumed to have the same temperature T+. A Bample 
multi -species ion curve with theoretical fit and results is shown in Figure 2. 


In the application of formula (1) to analyze ion I-V curves taken by a planar 
RPA, several approximations or assumptions are mode: 

(1) The ambient ionic velocity distribution is Moxwell-Boltzmann, i.e, , 
characterized by a temperature T+; and that this distribution is un- 
changed at the ItPA aperture, except for a uniform translation due to 
the satellite potential. 

(2) Each grid represents an equipotential plane according to the voltage 
applied to that grid, 

(3) The percentage of ions which pass through each grid is equal to the 
optical transparency of that grid taken at normal incidence. 

(4) Every ion which passes through all grids will hit the collector plate. 

(5) The measured current is due only to ions. 

In general, these approximations are not perfectly true. Consequently, the 
RPA should be designed and operated so as to maximize the validity of (1) - (5) , 
or else appropriate corrections should be made for effects which cannot be 
"designed away. ” 

Within the past few years, the operation of the RPA has been looked at in 
some detail, in order to better our understanding of its capabilities. These 
investigations have been undertaken for several reasons. For one, satellite- 
borne RPA's have been operating in altitude ranges where several ions (mainly 


0 + , H + and He 4 ) must be accounted for simultaneously. Also, the large Amount 
of data returned by a atoll' '.as and ground-based research stations has given a 
resonably good picture of the Ionosphere, so that experimenters are now con- 
cerned with accuracies of some percent, rather than orders of magnitude. 

Multi -experiment satellites have allowed intercomparison of results between 
different experimental methods. In addition, there have been disagreements 
between satellite experiments and radar backscattcr experimenters over iono- 
spheric temperatures (Hanson et al, , 1969; McClure and Troy, 1971), 

Some of the previously listed approximations under which RPA analysis 
occurs have been checked by different investigators: Whipple and Parker (1970) 
have considered the first approximation, and concluded that the ion I-V curve 
is minimally affected when the Debye length is small. The second approximation, 
that the grids are perfect equipotential planes, has been investigated by Knudsen 
(1966), Hanson et al., (1972), Golden et al. , (1972) and Whipple et al. , (1972). 
The last three assumptions were discussed by Whipple et al. , (1968), but only 
relative to high energy isotropic particles which pass through the RPA grids un- 
affected by the grid voltages, i.e. , not for the I-V curve. RPA theory for elec- 
trons has also been dealt with by Whipple and Parker (1969 a & b). 

It is the purpose of this paper to investigate the effects of grid transparency 
and finite collector size (approximations 3 and 4) on the determination of thermal 
ion density and temperature from I-V curves. The procedure will be to compute 


the shape of the I-V curve for given conditions (RPA dimensions, grid configura- 
tion, ion temperature and density, etc,), and then analyze the curve by the 
existing method. The differences between the assumed temperature-density 
values and the analyzed values will represent the magnitude of the error intro- 
duced by the grid transparency and finite collector effects. These results will 
not include the effect of grid plane potential non-uniformities (Whipple et al. , 
1972), 

FINITE COLLECTOR EFFECT 

Figure 3 illustrates the geometry used for determining the current to the 
collector. An ion of mass nij and velocity v arrives at the center of the aperture 
at polar angle 0 and azimuthal angle <f> (see figure 3a). The ion passes through 
the grids (not shown) and arrives at point P in the plane of the collector, P may 
be either on or off the collector. The ion generally does not follow a straight 
path from aperture to collector, due to the effect of electric fields between the 
grids. However, all positive ions with the same v, 0 , $ will follow similar 
paths, so that their locus in the collector plane will be a projection of the 
aperture, but displaced horizontally by a distance 1. The projection will lie 
either inside the collector (1 < R-r) so that all v, 9 , <p ions are collected; out- 
side the collector (1 > R+r) with no ions collected; or athwart the collector 
boundary (R+r > 1 > R-r) with some fraction f being collected. For a given 1, 
f is determined by 

f „ A (R, r, 1) 
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where A(R, r, 1) is the tire a of the projection within the collector boundary, 
i. e. , the intersection of the two circles . 

Since R and r are predetermined constants, f may be computed if 1 is known. 
Figure 3b illustrates the method of determining 1. The grids and collector are 
represented by potential planes, with a resulting electric field E between the 

mJ* 

planes. Between any two planes, E is assumed to be constant in magnitude and 
direction: 


where n is a unit vector normal to the potential planes, 

An ion enters this field at an angle 0j_j and velocity Vj_, . Its velocity 
parallel to the planes is unaffected, but its velocity normal to the planes changes 
depending upon E. The ion then passes plane j at an angle 0j ^ 0j„, and velocity 
, having traversed a horizontal distance 1 j (see fig. 3b). The values of 
0j, Vj and lj can be computed from 0j_j and Vj_j by simple mechanics. This 
calculation is done sequentially for all intervals d^ between the aperture and the 
collector, each ion entering the aperture with its ambient values of 0 and v. 

The value of 1 is then determined by summing 1 ^ for all intervals: 

1 = 21 . 
j 1 


For a given Ion mass and a given RPA configuration, lwill depend onfl and v, but 


will be independent of <f>. Therefore, the fraction f will be a function of 0 and v: 

f o f(0,v). 

The physical effect of tho finite collector size will bo dealt with in the cal- 
culation of I by using an effective aperture area of nr 2 f (0, v) , instead of the 
constant area wr 2 , 

TRANSPARENCY EFFECT 

The transparency of a square wire mesh grid has been derived by Whipple 
et al. , (1968): 
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where b is the wire diameter, a is the wire spacing and 0 , <t> are the ion angles 
at the grid (Fig. 3a). Since r (0, <p) is not strongly dependent on <f>, one can 
determine an average transparency (0) by integrating over <p : 

7T 

r(0) = “f 2 r(0,0)d0. 

7 TJ 
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The full expression is given by Whipple et al. , (1968). 

The transparency at each grid plane j is then given by r (0 j) , instead of the 
normal transparency r = r(0) used in equation (1) , The value of 0j depends 
upon the ambient 0, v of the ion before it passes the first grid. Therefore, we 
relable r (0j) as rj (0, v). Hie total transparency for an n-giid trap is then 

n n Ti (0, v) rather than r n (0). 
j = l J 
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CALCULATION OF THE ION I-V CURVE 


Wo can now incorporate the area factor f{0, v) and the transparency rj(0,v) 
into our calculation of the I-V curve. The standard I-V curve, equation (1), is 
calculated by integration ovor tho velocity distodbution function: 


N + -oo oo -1 |c x -v tx )2 +(v y -v |y ) 2 + (v z -v sz ) 2 l 

s 7ir 2 cr n (0) — — I dv x / dv v / dv. v, e !l L 


where v ml|1 is determined by tho retarding potential, and v sx , v sy , v JZ are the 
components of v in the RPA frame of reference. Changing to polar coordinates, 


we have: 
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the angles in polar coordinates for the velocity vector v s . Now we replace the 

area 7 rr 2 by tt r 2 f(0, v), and r 11 (0) by Il n r, (0, v). Doing this, and rearranging 

j « l J 


terms to the left, we have: 
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Formula (2) is now integrated numerically, and compared with the analyt- 
ically computed value of formula (1). For simplicity, formula (2) has been 
written for a single ion species. In practice, both formulas are normalized to the 
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ram ion current I 0 , which docs not depend on the retarding potential V: 

! 0 = rrr 2 T n (0)eN + v s . 

Naming the currents from formulas (1) and (2) as 1(1) and 1(2) respectively, we 
have: 


id) B F (s, t) S i+i e rf(s-t) + — ==e’< s -‘) 2 . 
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Both J(l) and J(2) depend on V: J(l) tlirough the term t 


A'(V-0 s ) 

- — , and 
k T\ 


J(2) tlirough the lower limit on the velocity integral shown In formula (2)* 
Therefore, formula (4) effectively describes the ion I-V curve including the 
finite collector and grid transparency effects, and formula (3) describes the 
I-V curve without the effects, i.e. , the planar approximation. 


RESULTS 

As a typical RPA, we use the one flown on ISIS-2, with r = 1,27 cm, R = 
3.17 cm, d « dj = d 2 - d 3 = 0,5 cm, and b/a = 0.04. To test the effect of 
finite collector size, runs ore made varying first R, and then varying d. The 
effect of increasing d will be similar to the effect of decreasing R: to cause 
more ions to miss the collector. However, an increase in d is not perfectly 
equivalent to a decrease in R except for r = 0, i.e., a point aperture. The 
results are presented as functions of normalized parameters d/r and R/r. Next, 
runs are made to determine the effect of grid transparency by varying b/a. To 
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separate tho effect of gild transparency from that of Unite collector else, wo 
sot the wire parameter b/a at os- s ^tinlly zero (ID* 4 ) while varying R/r and d/r, 
and we set R at essentially infinitely (10 motors) wliile varying b/a. All runs 
are made with the trap-velocity angle a a 0 . 

If all ions entered tho aperture at normal incidence, the planar approxi- 
mation F(s,t) would perfectly describe the shape of the current- voltage curve, 
insofar as the effects considered in this paper are concerned. But ions entering 
with high 0 and high v will see grid transparency different from r (0) , and will 
also have a chance of missing the collector. Low ion mass tmd high tempera- 
ture "will enhance the high 0 and liigli v component of the velocity distribution. 
Therefore, we choose hydrogen at 4000 °K. to give a reasonable upper limit on 
the difference between F(s,t) and J(?'» for c v<M RPA, i.e. , an upper limit on 
the effects we are investigating. 

Tlie accuracy of our calculation is tested by making a run with an RPA 
having b/a = 10 " 4 and R = 10 meters, where it is expected that J(2) “ F(s,t). In 
fact, the difference between the computed J(2) and F(s,t) for this case turns out 
to be less than 0.25%, which is to say that they are identical within the experi- 
mental accuracy of the RPA, 

The retarding voltage array uBed in the calculation is one point every 0..1 
volts from 0.0 to 1. 7 volts, which is enough to decrease the current by an order 
of magnitude. This amount of decrease is about what one might realistically 
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expect to observe for hydrogen experimentally; therefore ull our hydrogen 
calculations are based on a retarding voltage range of 0.0-1. 7 volts. 

Finite collector effect. 

Figure 4 shows the variation of the computed ion temperature and density 
with the collector parameter R/r. Other parameters ore fixed at b/a » 10" 4 
and d/r =0.39. The density Is normalized to the value obtained from the test 
trap (b/a = 10' 4 , R/r « 10 meters/1.27 cm = 787), and the assumed temperature 
is 4000°. As expected, the derived density decreases with decreasing R/r, so 
that there is a 30% density error when the collector is the same size as the 
aperture at R/r = 1. There is negligible error at the real RPA value of R/r = 
3.17cm/l. 27 cm = 2.5. 

The temperature, however, shows no truly significant error anywhere in 
the range Investigated, with a maximum error of about 1% (the error bars are 
assigned by the computer curve-fitting routine). We are left with the interesting 
result that the RPA can give the correct temperature even when it is giving 
density with a substantial error. 

The effect of varying is shown in Figure 5, for a trap with b/a = 10 -4 
and R/r = 2.50. The results bear out what was suggested earlier, namely that 
an increase in d/r is similar to a decrease in R/r, Again, the density error 
can be large (d/r > 1.0) while the temperature error is rather small . The real 
RPA value of d/r = 0.39 yields small errors in both density and temperature. 


'v’ 
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Grid transparency effect. 


Now that we have shown that the dimensions of the sample RPA are properly 
chosen, we investigate the effect of varying the grid transparency by varying 
b/a. Figure G shows these results for the RPA with a large collector (R/r » 
787). As with the finite collector effect in the previous section, the tempera- 
ture error remains small even for moderate errors in the density. The norm- 
alized density n + » however, shows a different dependency: n+ decreases linearly 
with b/a = 0, whereas n + approached 1.0 asymptotically as a function of d/r and 
R/r (see figures 4 and 0), The density error is therefore definitely non-zero, 
although it is less than 10% for most ion traps (b/a = 0.04 for the sample RPA, 
b/a< 0.1 for most other RPA's), 

It should be noted that the errors shown in the figures represent reasonable 
upper limits, being derived for hydrogen at 4000° K. Therefore, we check the 
lower limit by using oxygen at 1000° K for b/a = 0.2. The analysis results of 
this run (n+ shown in Fig. G, T + = 1014 +13° K) indicate that errors in the 
derived results for oxygen at 1000° K will be negligible over a reasonable 
range of b/a. 

SUMMARY 

We have investigated the effects of finite collector size and grid trans- 
parency on RPA ion data analysis under the assumptions outlined in the text. 
These effects result from the fact that ions generally do not enter the aperture 
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at normal incidence. The ions may lilt the grid wires with an increased proba- 
bility due to the variation of grid transparency with angle of incidence, or they 
may miss the collector entirely because of its finite size. Each effect decreases 
the current to the collector compared to that which is assumed by the planar 
approximation. 

Several current voltage curves are calculated for a varying UPA parameters 
and a given ion mass, temperature and density. The curves are analyzed under 
the planar approximation and the resulting densities and temperatures and 
plotted versus the varying parameter. The differences between the given and 
derived tomperatures and densities represent the errors due to the effects being 
investigated. 

We find that for the dimensions of our sample ltPA, the temperature and 
density errors are both small. The temperature error will remain small, even 
when the parameter is changed enough to substantially increase the density 
error. This information could be useful for someone designing an RPA to be 
used under conditions quite different from those of satellites in earth orbits. 

For instance, the aperture could be greatly increased for a given collector size 
to make measurements in low density regions. Or the grid wire size could be 
increased to combat heating effects during planetary probe missions. The de- 
rived temperature can be assumed correct, and the density increased according 

« 

to calculations. 
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FIGURE CAPTIONS 


Figure 1. 

Figure 2. 
Figure 3a. 

Figure 3b, 
Figure 4, 

Figure 5. 

Figure 6. 


Schematic of the RPA, showing the aperture radius r, the collector 
radius R, and the trap depth dsdj+dj + dg. 

A multi-component current-voltage ourve, with analysis results. 

The projection of the aperture image onto the collector plane by 
ions having specific velocity magnitude and direction, 

» 

Variation of the ion path direction between two grids. 

Variation with collector size of ion temperature and density derived 
from computed I-V curves for an RPA with d/r =0,39 and b/a = 
10" 4 , All points are for H + at 4000° K. 

Variation with RPA depth of ion temperature and density derived 
from computed I-V curves for an RPA with R/r = 2,50 and b/a = 
10 -4 . All points are for H* at 4000° K. 

Variation with wire parameter b/a of ion temperature and density 
derived from computed I-V curves for an RPA with R/r = 787 and 
d/r = 0.39. The triangle is for 0 + at 1000 °K. ; all other points 
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